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ABSTRACT 

We investigate the formation and evolution of nested bar systems in disk galaxies in a cosmological 
setting by following the development of an isolated dark matter (DM) and baryon density perturbation. 
The disks form within the assembling triaxial DM halos and the feedback from the stellar evolution is 
accounted for in terms of supernovae and OB stellar winds. Focusing on a representative model, we 
show the formation of an oval disk and of a first generation of nested bars with characteristic sub-kpc 
and a few kpc sizes. The system evolves through successive dynamical couplings and decouplings, 
forcing the gas inwards and settles in a state of resonant coupling. The inflow rate can support a 
broad range of activity within the central kpc, from quasar- to Seyfert- types, supplemented by a 
vigorous star formation as a by-product. The initial bar formation is triggered in response to the tidal 
torques from the triaxial DM halo, which acts as a finite perturbation. This first generation of bars 
does not survive for more than 4-5 Gyr: by that time the secondary bar has totally dissolved, while 
the primary one has very substantially weakened, reduced to a fat oval. This evolution is largely due 
to chaos introduced by the interaction of the multiple non-axisymmetric components. 
Subject headings: galaxies: evolution - galaxies: formation - galaxies: halos - galaxies: kinematics 
and dynamics - galaxies: starburst - galaxies: structure 



1. INTRODUCTION 

The formation and evolution of galaxies are closely as- 
sociated with the radial redistribution of baryonic and 
dark matter (DM) and of their angular momentum (J) 
(e.g., Athanassoula 2002, 2003). This process is further 
amplified by the lack of an axial symmetry in the ba- 
sic galactic components over substantial periods of time. 
Numerical simulations have shown that DM halos form 
universally triaxial, i.e., flattened and elliptical (e.g., 
Allgood et al. 2006), and remain in this state for the 
few Gyr during the disk growth (Berentzen & Shlos- 
man 2006). A clear majority of disks are barred in the 
nearby universe (e.g., Knapen et al. 2000; Grosbol et 
al. 2004; Marinova & Jogee 2006), and a fraction of 
barred disks is maintained at least to z ~ 1 (Jogee et al. 
2004; Elmegreen et al. 2004; Sheth et al. 2003). Simula- 
tions have also shown that galactic disks are subject to 
the spontaneous and induced stellar bar formation (e.g., 
Athanassoula 1984). A wide range of studies, theoreti- 
cal and observational, argue that disk galaxies spend a 
substantial fraction of their life in a non-axisymmetric 
stage. Their morphological components are subjected to 
mutual gravitational torques that, in addition to exter- 
nal factors, provide an efficient mechanism for driving 
their internal evolution (e.g., Weinberg 1985; Athanas- 
soula 1992; Heller & Shlosman 1994; Sellwood 2006). 

The shapes of the host halos affect strongly the grow- 
ing disk, especially because early disks are dominated 
by the dissipative baryonic component, unlike disks in 
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the nearby universe. The gas responds dramatically to a 
non-axisymmetric driving by shocking and loses its rota- 
tional support. One can expect that bars or grand-design 
spiral arms will be triggered in the disk and that the 
central mass concentration will build up. The most effi- 
cient redistribution of J, down to the smallest radii, in a 
disk can be attributed to the mechanism of nested bars, 
stellar and gaseous, the so-called bars-in-bars scenario 
(Shlosman et al. 1989, 1990). The system of nested bars 
that consists of the primary (large-scale) and secondary 
(sub-kpc) bars tumbling with different pattern speeds, 
can facilitate the radial inflow of gas and fuel the ac- 
cretion processes onto supermassive black holes (SBHs) . 
The latter correlate tightly with the properties of galac- 
tic bulges across seven decades in radius at least (e.g., 
Ferrarese & Ford 2005). About 1/3 of barred galaxies 
have secondary bars (Laine et al. 2002; Erwin & Sparke 
2002). Furthermore, formation of the SBH can be trig- 
gered by this process, because it leads to a specific en- 
tropy minimum in the center (Begelman et al. 2006). 
The dynamics of nested bars has been investigated (e.g., 
Shlosman et al. 1989; Friedli & Martinet 1993; Combes 
1994; Maciejewski & Sparke 2000; Heller et al. 2001; 
Shlosman & Heller 2002; El-Zant & Shlosman 2003), but 
modeling of nested bars in a cosmological scenario even 
for isolated DM halos was never attempted or validated. 

This Letter demonstrates how a nested bar system 
forms in a disk growing in an assembling triaxial DM 
halo. We use cosmological initial conditions and follow 
the Hubble expansion and the subsequent collapse of an 
isolated perturbation in the gas and the DM. Our star 
formation (SF) algorithm is physically motivated and we 
include feedback from stellar evolution. A large num- 
ber of models have been run and show the formation of 
galactic disks whose structural parameters fit within the 
observed range. Here we describe a representative model 
only and focus on the formation of nested bars. Addi- 
tional aspects will be addressed elsewhere. 
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2. NUMERICS, STAR FORMATION, FEEDBACK AND 
INITIAL CONDITIONS 

The simulations were performed with an updated ver- 
sion of the FTM-4.4 hybrid A-body/SPH code (Heller & 
Shlosman 1994), using the routine f alcON (Dehnen 2002) 
to compute the forces. The gas temperature was ob- 
tained from the energy equation. We used Nn M = 5 x 10 5 
and -/Vsph = 5 x 10 4 , gravitational softening of 150 pc 
for DM and stars, and dynamic softening with a mini- 
mum of 250 pc for the gas. Tests have been performed to 
check the sensitivity of the results to the algorithm and 
its parameters. A pure DM run conserved the energy to 
within 1% and J to within 0.1%. 

The modified prescription for SF and the feedback from 
stellar evolution take place in the Jeans unstable, con- 
tracting region. We fix the gas-to-background density at 
which gas is converted to a star, the local collapse-to- 
frcc fall time, and introduce the probability that the gas 
particle produces a star during a given timestep. Four 
generations of stars form per gas particle with an instan- 
taneous recycling along with an increment in metallicity. 
The balance of the specific internal energy along with 
the gas ionization fractions of H and He and the mean 
molecular weight are computed as a function of p and T 
for an optically thin primordial composition gas. Feed- 
back from stellar evolution includes the supernovae and 
OB stellar winds, and uses the thermalization parameter 
— a fraction of the feedback deposited as a thermal en- 
ergy, and converted into kinetic energy via equations of 
motion (Heller et al., in prep.). 

The initial conditions are those of a spherically- 
symmetric density enhancement at z — 36 evolved using 
an open CDM (OCDM) model with fl a = 0.3, h = 0.7. 
The difference between the OCDM and ACDM is min- 
imal on sub-galactic scales. The initial density profile 
corresponds to the average density around a 2a peak in 
a Gaussian random density field. The spin parameter 
A = 0.05 is set by the angular velocity u> oc r _1 , where 
r is the cylindrical radius and the central kpc has been 
softened. We impose the collapse redshift of z = 2. The 
mass of collisionless DM particles is 7 x 10 11 M and the 
gas comprises 10% of the total mass initially. 

The main simplification in our models is that we ne- 
glect interactions and minor or major mergers. Yet, be- 
cause the evolution is determined to a large extent by the 
nonlinear interaction between three non-axisymmetric 
components, the results presented here can be consid- 
ered representative. 

3. RESULTS 

The initial perturbation expands with the Hubble flow 
and collapses, forming a DM halo — a nearly nonrotating 
triaxial figure. In particular, the innermost 10 kpc ex- 
perience a very abrupt and short-lived increase of their 
ellipticity to e p 0.2, presumably due to radial orbit 
instability. The disk grows within the halo equatorial 
plane, visible already at r ~ 0.5 Gyr. Its growth nearly 
washes out the inner, < 10 kpc, halo ellipticity, due to 
the increased central mass concentration in the system 
and the out-of-phase response by the disk. The inner 
5 kpc (10 kpc) of the DM become nearly axisymmetric 
and only slightly flattened after ~ 2 Gyr (~ 7 Gyr), while 
at larger radii e p — » 0.15 — 0.25 and f p — > 0.3 — 0.4. 

The disk evolution consists of gas-dominated (first 




Fig. 1. — Dynamical evolution of nested bars — a snapshot at 
r ~ 3 Gyr (see also the Animation Sequence 1). A large-scale bar 
(left frame) is shown in gas (green), stars (yellow) and SF (blue), 
and the nuclear bar in the same colors (right frame) . The left frame 
is 24 kpc and the right frame is 4.8 kpc on the side. The Animation 
shows the evolution of the disk and the large-scale primary bar in 
the above colors (left frame) and of the inner kpc and the associated 
secondary nuclear bar (right frame) over 4 Gyr. 



1.5 — 2 Gyr) and star-dominated phases (Fig. 1 and An- 
imation Seq. 1). The ratio of baryonic-to-DM matter 
within the central 10 kpc increases from ~ 0.4 to ~ 0.54 
over the Hubble time due to an adiabatic contraction. 
The corresponding pure DM model results in the NFW 
(Navarro et al. 1996) density profile and the character- 
istic radius i? n fw ~ 9 kpc and is even more triaxial. For 
the model with baryons — R n fw ~ 4.5 kpc, and the fit 
quality has worsened. The cusp is baryon-dominated. 

Initially, the disk is roughly oval and dominated in the 
outer parts by m = 2 and 3 grand-design spirals. Grad- 
ually it becomes more axisymmetric. The characteristic 
time for the disk (and halo) buildup, i.e., reaching 50% 
of its mass at z = 0, is ~ 1.5 Gyr. The rate of a baryonic 
inflow into the disk region crests within 1 Gyr. The SF 
rate shows a similar behavior. 

The disk grows from inside out. A primary ~ 2 kpc- 
radius bar appears in the stellar component around 
0.6 Gyr normalto the halo major axis, increasing to 3 kpc 
by 0.7 Gyr. By 0.8 Gyr, the inner part of the 4 kpc bar, 
which is gas-dominated, decouples in the prograde direc- 
tion and forms a secondary bar. By 1.5 Gyr, most of the 
gas in the secondary has fallen to the center — further 
collapse is inhibited by the numerical resolution. 

Fig. 2a follows the evolution after 1 Gyr. Both bars 
are strong at this time, as given by A2 ~ 0.6 and 0.45 — 
the amplitudes of corresponding m — 2 modes. They 
are SF-dominated and involve a large fraction of the 
stellar and gas disk mass. The shape of the primary 
depends strongly on its orientation with respect to the 
halo and disk major axes — main sources of tidal torques 
— and on its internal (self-gravitational) response, e.g., 
between 1 to 2 Gyr (Anim. Seq.). The disk response 
to the halo torques drives a pair of strong grand-design 
arms whose pitch angle gradually decreases forming a 
pseudo-ring outside the bar, after which the open arms 
are regenerated (Anim. Seq.). The bars appear to decay 
after ~ 1.7 Gyr in tandem with the arms, but strengthen 
abruptly after 2 Gyr. The primary weakens again after 
2.5 Gyr, while the secondary bar dissolves by ~ 5 Gyr. 
The bar sizes vary between r s ~ 0.5 — 1 kpc for the sec- 
ondary and r p ~ 5 — 7 kpc for primary bars. 

The central issue of bars-in-bars systems is the dynam- 
ical coupling of bars as measured by the ratio of pattern 



Induced Nested Bars in Assembling Halos 



3 




Fig. 2. — Nested bars evolution: (a) m = 2 amplitudes A2, and (b) associated pattern speeds. Shown are the stellar A2 and pattern 
speeds of the secondary (red line) and primary (green line) bars and of the oval disk, between 5 — 10 kpc (blue line). 
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Fig. 3. — Gas inflow rate across the central kpc. 



speeds, fl s /fl p (Fig. 2b). The primary bar slows down 
prior to 1.5 Gyr, tumbles with a constant f2 p over the 
next 1 Gyr, speeds up and stabilizes around 2.8 Gyr 
for the rest of the simulation time. The secondary bar 
speeds up dramatically before 1.5 Gyr, stabilizes for 
~ 0.2 Gyr, then slows down and stabilizes again after 
2.8 Gyr. Its fi s appears to have a plateau also around 
2.4 Gyr. Hence the bars go over a number of stages with 
various f2 s /f2 p « const. 

Additional evidence for coupled evolution comes from 
varying corotation (CR) radii, r cr s and r crjP , of bars and 
from the position of the inner Lindblad resonance (ILR) 
of the primary, nLR, p . During the first 1.5 Gyr, r crjS 
and r crjP increase sharply from ~ 1 kpc and 4.2 kpc to 
2.3 kpc and 8.2 kpc, respectively. This increase is caused 
by a substantial inflow that leads to a surge of the rota- 
tion curve. This is followed, between 1.5 — 2.8 Gyr, by a 
plateau in r cr . s (T) and r CI%p (r), and a further increase of 
r cr , s to 3.5 kpc. Both remain stable thereafter. We esti- 
mate r cliP /r p ~ 1.1 — 1.6 and r cliS /r s ~ 1.3 — 5. The value 
of 5 is achieved at the end of the first 1.5 Gyr (it results 
from the runaway action of the gas-dominated secondary, 
see §4), and after 2.8 Gyr — making the secondary much 
shorter than its CR. The ILR of the primary can be in- 
ferred to lie outside the nuclear ring (see below) during 
this stage, riLR iP ~ 3.5 kpc. 

Nuclear rings form at the interface between the bars, 
and a pair of grand-design arms extend along the pri- 
mary — the ring and the arms are delineated by the 
SF intermittently. The ring evolution is more dramatic 
when viewed in the SF colors and it fades away with the 
weakening of the primary. 



4. DISCUSSION 

We have simulated, in a large number of models, the 
formation of galactic disks in assembling triaxial DM ha- 
los and analyzed a representative model. The accompa- 
nying SF and the stellar feedback have been included. 
We followed the collapse of an isolated cosmological den- 
sity perturbation with A = 0.05 from its linear regime, in 
the OCDM universe. The disk and halo form, i.e., reach 
50% of their final mass, over a period of ~ 1.5 Gyr. The 
DM halo develops a clear triaxial shape, elliptical in the 
plane perpendicular to the original angular momentum 
axis and flattened along this axis. The resulting halo 
figure tumbles very slowly, ~ n over the Hubble time. 

The halo triaxiality decreases during the disk growth 
because of the increasing central mass concentration and 
the disk response to the halo potential. The first effect 
reduces somewhat f p and e p of the halo, even in the ab- 
sence of a baryonic component. The second effect results 
from the negligible tumbling of the halo figure — the 
inner ILR and the outer ILR are pushed to the center 
and to large radii, respectively. The disk responds out- 
of-phase with the halo potential, thus diluting the halo 
ellipticity in its equatorial plane. This effect is not re- 
lated to dissipation. 

Most importantly, the halo torques lead to a strongly 
oval-shaped, growing disk and trigger the bar forma- 
tion — reminiscent of torques exerted on the disk during 
galaxy interactions. These initial bars form as a result of 
a finite amplitude perturbation and not as a consequence 
of an exponential growth from an infinitesimal perturba- 
tion. Their formation timescale is much shorter than for 
spontaneous bars and they appear first normal to the 
halo major axis. In our simulation the initial triaxial- 
ity of the inner halo is due to a short-lived episode of 
the radial orbit instability erasing the initial conditions 
(e.g., MacMillan et al. 2006). Other causes of triaxiality 
are possible — we wish to stress here its effect on bar 
formation and evolution. 

While nested bars have been simulated before, with 
various degrees of self-consistency, this simulation differs 
in at least four major aspects. First, this is the most 
self-consistent model of nested bars in the literature — 
no a priori assumptions have been made. Second, this 
is the first model where nested bars form from cosmo- 
logical initial conditions. Third, both bars are very gas 
rich, while previous models of nested bars focused on 
purely stellar or gaseous bars (Shlosman 2005, for re- 
view) . Fourth, we have incorporated SF and feedback — 
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the stellar and gaseous masses in the simulations are not 
individually conserved. 

It is helpful to divide the nested bar evolution into 
three distinct phases, namely formation, dynamical ad- 
justment and dynamically quiescent phases, as exhib- 
ited in Fig. 2b. The initial phase, which is preceded 
by an early disk, ends at ~ 1.5 Gyr and extends over 
~ 0.7 Gyr from the time the secondary can be identified. 
The primary and especially secondary bars form gas rich. 
During this time, f2 p decays steadily while f2 s increases 
monotonically by a factor of ~ 2. The latter increase 
is due to the gas collapse by a factor of 3-4 within the 
secondary, down to numerical resolution scales. Simula- 
tions of a dynamical runaway of a pure gaseous secondary 
have shown f2 s ~ r^T 1 (Englmaier & Shlosman 2004). A 
smaller growth in Cl s in our model can be explained if a 
correction is made for the non-dissipative stellar compo- 
nent in the bar. Note that f2 p follows closely the pattern 
speed of the oval disk, which forms in response to the 
halo ellipticity (Fig. 2b). 

The second phase of evolution, between 1.5 - 2.8 Gyr, 
is characterized by a dynamical adjustment of primary 
and secondary bars via their pattern speeds. The ratio 
f2 s /f2 p ~ 6.1 is steady for about 0.2 Gyr, then decays 
steadily because of a monotonic decrease in f2 s , while f2 p 
stays constant initially, then increases slowly. The bars 
pass through an intermediate period of ~ 0.2 Gyr when 
their pattern speeds are possibly locked, Q s /Q p ~ 4.4, at 
about 2.4 Gyr. A dramatic inflow leads to a buildup of 
a massive gas disk that generates grand-design shocks in 
response to the (elliptical) halo driving when the disk el- 
lipticity surpasses a critical value (B. Pichardo & I. Shlos- 
man, 2005, unpublished). This process largely deter- 
mines ftp which is also the pattern speed of the spiral 
arms. Hence, contrary to other models, the overall evo- 
lution in the system is driven by the halo triaxiality. 

While the disk shape depends strongly on the angle 
with the halo major axis, the primary bar is shaped by 
the angle with the halo and with the oval disk — an 
efficient way to amplify chaos within the bar, weaken- 
ing it on a short timescale (El-Zant & Shlosman 2002; 
Berentzen et al. 2006; Berentzen & Shlosman 2006). 
This can explain the sharp decay of primary amplitude 
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at 1.7 Gyr and 2.5 Gyr, both accompanied by massive 
spiral arms in the disk. Secondary bars have been found 
to adjust their properties, such as axial ratios and radial 
extent, depending on the mutual bar orientation (Heller 
et al. 2001; Shlosman & Heller 2002). The situation is 
much more complex in the present run, because a num- 
ber of asymmetric components coexist. The individual 
interactions, which are not always possible to disentan- 
gle, will increase the fraction of chaotic orbits in the bars 
and generate local (non-grand-design) shocks in the gas. 
These interdependencies can have observational corollar- 
ies at redshifts corresponding to disk growth and will be 
addressed elsewhere. 

The third phase, after ~ 2.8 Gyr when the gas fraction 
is low, is that of a stable coupling between fl p and f2 S) 
i.e., f2 s /f2 p ~ 2.7, confirming that stellar nested bars 
can be long-lived systems (El-Zant & Shlosman 2003). 
Because hlr.p ~ ?*cr,s at this time (see §3), it is possible 
that we observe a resonant coupling between primary and 
secondary as proposed by Tagger et al. (1987) for any 
two modes. 

During the first stage, the gas inflow across the cen- 
tral kpc, Mi, corresponds to a quasar-type activity of 
60 M yr" 1 (Fig. 3), while the SF peaks at 25 M yr" 1 . 
In the second stage, Mi crests at ~ 25 M Q yr _1 around 
2 Gyr, dropping to its lowest value close to zero there- 
after. Nuclear rings and the associated SF are intimately 
related to this inflow (Athanassoula 1992; Heller & Shlos- 
man 1996; Knapen 2005), but in nested bars they vary 
in shape and have a more limited life span in response to 
the time-dependent potential (Shlosman & Heller 2002). 

Hence, we have demonstrated that nested bars in iso- 
lated halos form from cosmological initial conditions and 
go through a series of dynamical couplings and decou- 
plings, while channeling the gas inwards to the small- 
est scales resolved numerically. This inflow can support 
the early quasar-type and Seyfert-type activity there- 
after. The exact conditions leading to an SBH formation 
(Begelman et al. 2006) are beyond the scope of this work, 
but the remote 'boundary' conditions are in agreement. 
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